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We have used a study of the process e'^e~ — » fj,'^ ^~"f at a center-of-mass energy near the T(4S) 
resonance for a fj,'^ j-L~ invariant mass range near the J/ijj mass to extract the cross section a{e^e~ 
J/^l ~* t)- The data set, corresponding to an integrated luminosity of 88.4 fb^^, was 

collected using the BABAR detector at the PEP-II collider. We measure the product r(J/-0 — > e'^e^)- 
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B{J/ip m^M") to be 0.330±0.008±0.007 keV. Using the world averages for B{J/ip and 
B{J/ip — > e^e~), we derive the J/ip electronic and total widths: r(J/i/) e^e~) = 5.61 ±0.20 keV 
and r = 94.7 ± 4.4 keV. 

PACS numbers: 13.20.Gd, 13.66.Bc, 14.40.Gx 



The possibility of using e~^e~ annihilation with ini- 
tial state radiation (ISR), e+e^ ^ hadrons + 7, to mea- 
sure the e^e^ cross sections into hadrons over a wide 
range of center-of-mass (CM) energies in a single ex- 
periment has been discussed in the literature [jj]. In 
this paper, we have implemented this idea by study- 
ing the process e~^e~ — > fj,'^^~"f for fi^fJ.~ masses in 
the range from 2.8 to 3.4 GeV/c^. We measure the 
cross section <T(e+e^ J/il'l ^ /^^M~7) ^'^'^ derive 
the product of electronic width times branching fraction 
T{J/i]] e+e~) • ^(J/t/i M"*"/^")- Using the world 
averages for B{J/tp fi^fi^) and B{J/tp e+e~) 0, 
we then derive the electronic and total widths of the J/ip 
meson. The data used in this analysis were collected with 
the BABAR detector at the PEP-II asymmetric e+e" 
storage ring 0. 

The Born cross section for the process e+e^ —>■ fi'^ fi^j 
in the J/ip mass region has contributions from three 
Fcynman diagrams, as illustrated in Fig.^ The first and 
second of these diagrams describe the pure QED pro- 
cesses corresponding to initial state radiation and final 
state radiation (FSR). The visible QED cross section in 
BABAR (defined by our ISR photon acceptance) is about 
1.2 pb in the di-muon mass range 2.8-3.4 GeV/c^. The 
contribution of the FSR process to the QED cross section 
depends on the photon energy and angle, and is about 
10-20% for the kinematic regime we study. The interfer- 
ence between ISR and FSR amplitudes does not change 
the total cross section, but leads to charge asymmetries 
in the muon angular distributions. 

The Born cross section for J/i/' production (Fig.^) is 
given by 



da; 



W{s,x)-ao{s{l-x)), 



s, 



(1) 



where ^/s is the e"'"e~ invariant mass, x = 2E^ 
is the photon energy in the CM, and ctq is the Born cross 
section for e+e~ J/tp fi^ fi^ . The function 



2a 



M^(s,x) = — •(21n^-l).(l-x+-) (2) 

TTX me Z 

describes the probability of ISR photon emission. Here 
a is the fine structure constant and Trig is the electron 
mass. ISR photons are emitted predominantly at small 
angles relative to the electron direction. About 10% of 
the photons have CM polar angles in the range 30° < 
9 < 150° and can be detected in BABAR. 

As a first approximation, the Born cross section for 
e^e^ — > J/ijj —^ is given by the Breit-Wigner for- 




FIG. 1: Diagrams for e^e 7- a) Initial state radia- 

tion, b) Final state radiation, c) J /ijj production. 
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where m and F are the J/ip mass and total width, 
respectively, B^e and i?^^ are the branching fractions 
B{J/'ip e^e~) and B{J/ip —* For a narrow 

resonance, such as the J/^p- we can replace the Breit- 
Wigner with a S function 7rmFi5(s — m?) and integrate 
over photon energy to find 
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Here Fee = T • B^e- These formulae do not account for 
interference between the e^e~ ^ J jip ^ A*"*"/^" and the 
non-resonant (QED) e+e~ amplitudes. 

The cross section for e^e~ — > /x"*"^" including QED 
and resonant J/i/' production amplitudes and their in- 
terference can be written as 
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mT 



imT 



(5) 



where Q ~ 'i^jB^^^W^/ a. The interference term changes 
sign at the J/ip mass. Therefore, it does not change the 
integrated cross section of Eq. significantly, but does 
change the shape of the mass distribution. Because the 
J /ip cross section is so much greater than the QED cross 
section at resonance (Q^ « 600), the power-law behavior 
of the Breit-Wigner tails produces observable interference 
even 1000 widths from resonance. The expected di-muon 
mass spectrum, convolved with the detector resolution, is 
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FIG. 2: The di-muon mass spectrum calculated with (solid 
line) and without (dashed line) interference between the reso- 
nant J/tl} production and QED amplitudes after convolution 
with the detector resolution function. 

shown in Fig. |21 The interference is clearly seen, despite 
the experimental resolution, 14.5 MeV/c^, being more 
than 100 times the J/ip natural line width. The max- 
imum relative difference between the spectra calculated 
with and without interference is about 7%. The interfer- 
ence also leads to a 1.3 McV/c^ shift between the max- 
imum of the resonance peak and the actual J/ip mass. 
The shape of the expected mass spectrum is very sensi- 
tive to the tails of the Breit-Wigner approximation used 
in Eqs. ^ and (0), where its validity far from resonance 
is questionable. To estimate the sensitivity of our analy- 
sis to the details of the shape assumptions, we will take 
the full difference between fits that do, and do not, use 
interference as a measure of the systematic uncertainty. 

The width of the J/iJj has been measured directly in 
pp annihilation with the result 99 ± 12 ± 6 keV 0. In 
e~^e~ annihilation, measuring the area under the reso- 
nance curve for e^e~ s- J/tJj fJ-^ fJ-^ gives the prod- 
uct Fee ■ Bfj^f^ as seen in Eq. Combining this with 
the leptonic branching ratio yields the total width. The 
BES Collaboration made a comprehensive collection of 
measurements at the J/i/' from which they determined 
F = 84.4 ± 8.9 keV 0. This superseded results ob- 
tained from original measurements made of the area un- 
der the excitation curve in 1975. More recently, the BES 
Collaboration has measured the leptonic branching ra- 
tio with a 1.5% uncertainty using J/i/j's from the decay 
ip{2S) J/TpTTTT 0. It is this result that we combine 
with our measurement of Fee ■ S^i/^ to obtain the highest 



precision result to date for the total width of the J/ip. 

Charged particle tracking for the BABAR detector is 
provided by a five-layer silicon vertex tracker (SVT) and 
a 40-layer drift chamber (DCH), operating in a 1.5 T ax- 
ial magnetic field. The transverse momentum resolution 
is 0.47% at 1 GeV/c. Energies of photons and electrons 
are measured by a CsI(Tl) electromagnetic calorimeter 
(EMC) with resolution of 3% at 1 GeV. Charged particle 
identification is provided by ionization measurements in 
the SVT and DCH, and by an internally refiecting ring- 
imaging Chercnkov detector (DIRC). Muons are iden- 
tified in the solenoid's instrumented flux return (IFR), 
which consists of iron plates interleaved with resistive 
plate chambers. The data sample used for this analy- 
sis corresponds to an integrated luminosity of 88.4 fb""'^ 
recorded in the vicinity of the T{AS) resonance. 

The initial selection of fj,'^ candidates requires that 
all particles are detected inside a fiducial volume and that 
the event kinematics are consistent with the hypothesis 
e~^e~ — *■ /i^/i~7. Photons must have polar angles in the 
range 0.35 < 6* < 2.4 radians and must have a CM energy 
above 3 GeV. Muon candidates must have polar angles in 
the range 0.35 < 9 < 2.4 radians and transverse momenta 
above 0.1 GeV/c, and must originate from the interaction 
point. Energy and momentum balance is provided by the 
conditions |i?totai - £^boams| < 1-5 GeV and A^P < 0.07. 
Here ii'totai is the summed energy of the muon candi- 
dates and the photon, A^P is the angle between the pho- 
ton and the direction of the di-muon missing momentum, 
Pniiss = Pe+ + Pe- ^ ^ P^- • Wc reduce backgrounds 
using a one-constraint fit to the hypothesis that the recoil 
mass against the di-muon be zero. Requiring < 20 re- 
jects 90% of the multihadron ISR contamination (general 
e+e~ — > qq^y reactions) and about 10% of signal events. 

The large background from e~^e~ — > e'^e~j is sup- 
pressed by requiring that the charged track momenta 
be greater than 0.5 GeV/c and that the corresponding 
energies detected in the calorimeter be small: E'emc < 
0.4 GeV for each track. The average energy deposition 
of muons in the calorimeter is about 0.2 GeV, while elec- 
trons typically deposit more than 90% of their energy 
in the calorimeter. Additional suppression of this back- 
ground is achieved by requiring large angular separation 
(in the CM) between the charged tracks and the pho- 
ton (cos6'*^ < 0.5). This also reduces the level of FSR 
/i"^/z~7 events in the final sample by a factor of two. The 
invariant mass distribution of approximately 70000 di- 
muon pairs in our final sample is shown in Fig. 01 About 
7800 events are in the J/ip peak. 

To increase our detection efficiency and minimize sys- 
tematic uncertainties, we do not use IFR information in 
selecting muons. However, we do use this information for 
estimating backgrounds. For this purpose, muon identi- 
fication requires that a track penetrate at least 2 nuclear 
interaction lengths (A) of IFR material, and that the dif- 
ference between the measured and expected muon ranges 
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be less than 2A. This algorithm is 90% efficient for true 
muons and misidcntifics about 10% of real electrons as 
muons. 

The remaining electron contamination in our final sam- 
ple is estimated using a subsample of events enriched with 
electrons. We require that neither muon be identified in 
the IFR using the algorithm described above. We then 
require that the DCH based dE/dx measurements for the 
two tracks be consistent with the di-electron hypothesis 
and inconsistent with the di-muon hypothesis. This elim- 
inates 95% of di-muons and retains 85% of di-electrons. 
The fraction of electron events in the final sample is es- 
timated to be (0.1 ± 0.1)%. 

ISR events with hadronic final states are another 
source of background, both on resonance and off. For 
e+e" J/i^l 7r'*"7r~7 and e~^e~ J ^ 
K^K^-j, the cross sections are proportional to the ra- 
tios of branching fractions B{J/ip tt^tt^)/B{J/iP 
H+ti-) w 2.5 • 10^3 and B(J/V' ^ K+K-)/B{J/iP 
jX^fx") K, 4 • 10^^. To first approximation, the off- 
resonance ratios, a{e^e'^ —>■ 7r+7r^)/(T(e+e^ —>■ ^^fi^) 
and a{e~^e~ K~) / a{e^ e" fi^ii~) are similar 

to those on resonance. As off-resonance production pro- 
ceeds via virtual photon intermediate states while on- 
resonance production proceeds via both virtual photon 
intermediate states and hadronic intermediate states 0, 
the on-resonance ratios are overestimates of the off- 
resonance ratios. Thus, we consider the on-resonance 
background rate as an upper limit for both. 

The suppression of kaon and pion reactions was stud- 
ied using samples of e+e^ ujj — + Sttj and e+e^ 
(jyy K~^K~'y events. About two thirds of these events 
are rejected by the calorimeter energy deposition require- 
ment. Under the di-muon hypothesis, the peak of the 
J/ip ^ K^K~ distribution transforms into a broad dis- 
tribution with A/p^ < 2.95 GeV/c^. The only back- 
ground peaking under the J/ip — > signal is that 
due to J/f/' — > 7r+7r~; its contribution is estimated to be 
(0.09 ± 0.03)%. The only other decay into two charged 
hadrons, J/ip pp, produces events with M^^ < 2.4 
GeV/c^, and thus contributes no background in the di- 
muon mass range studied here. The total non-resonant 
background from e^e^ —>■ e^e~^, 7r"'"7r~7, _ftr+iir~7 pro- 
cesses is estimated to be (0.3 ± 0.2)%. 

The background from ISR production of higher multi- 
plicity multihadron events is estimated from Monte Carlo 
simulation. We estimate the background from multi- 
hadron J/i/; decays to be less than 0.05% using simu- 
lated e+e^ J/ipJi J/'fp ~^ Stt events and J/ip charged 
particle multiplicity data j^. We also note that such 
events populate the mass region below 3 GeV/c^ when 
misidentified as signal events. We use the JETSET 
event generator to simulate the hadronic part of the 
e+e^ qqj, q = u,d,s cross section. We find the back- 
ground due to such events to be less than 0.3%. As these 
background rates are not the dominant sources of systcm- 
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FIG. 3: The mass spectrum for observed events. The curve 
is the result of the fit described in the text. 

atic uncertainty in our final results, we have not tried to 
determine them with greater precision. 

We use a binned maximum likelihood fit to describe 
the mass spectrum of Fig. 13 The mass range used, 2.8- 
3.4 GeV/c^, is divided into 150 bins of width 4 MeV/c^. 
The probability density function (PDF) for the J/ijj sig- 
nal is modeled as the convolution of a J/i/i Breit-Wigncr 
line shape and the resolution function shown in Fig. 0] 
This is derived from detector simulation in conjunction 
with an e~^e~ — *■ ^~^fj.~"f event gene rator based on the 
differential cross sections of Ref. [ll| . Soft photon radi- 
ation is generated with the use of the structure function 
method of Ref. [l^ and the PHOTOS package |l3 for 
electron and muon bremsstrahlung, respectively. Muon 
bremsstrahlung leads to the low mass tail observed in the 
spectrum of Fig. 01 To account for possible resolution 
differences between simulation and data, the resolution 
function shown in Fig. 0] is convolved with an additional 
Gaussian smearing function of width ac- Both ctg and 
the observed J/^ peak position, Mj/^p, are parameters in 
our fit. A Monte Carlo calculation shows that the shape 
of the non-resonant cross section can be described well 
by a linear function. To account for possible deviations 
from this hypothesis {e.g., due to detector response) a 
second-order polynomial is used to fit the experimental 
spectrum. The full PDF is written as 

No 
L- (nii ) 

1 + a{m, - Mj/^) + h(m, - Mj/^)% (6) 

where rrii is the central value of the zth bin of the data 
histogram, Nq = ^ • Am is the level of the non-resonant 
mass distribution at m = Mj/^, Am = 4 MeV/c^ is the 
bin width, H is the PDF for the J/ 4> signal with detector 
resolution, and a and h are the background polynomial 
coefficients. To account for the interference between res- 
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FIG. 4: The distribution of the reconstructed mass minus 
the generated mass in Monte Carlo events. 

onant and non- resonant amplitudes described in Eq. 
the PDF is divided by the correction function C{mi), 
which is the ratio of the di-muon mass spectra calcu- 
lated with and without interference as shown in Fig. |2] 
Because the shape of this function depends on the J/ip 
parameters (mass, full width), an iterative procedure is 
used to calculate it. The ratio 

(7) 



R 



•j/i, 



^ • Am 

dm 



is the main fit parameter. Here Nj/^ is the number of 
observed J/V' decays. After substituting cross sections 
for numbers of events, this ratio can be rewritten 



R 



'J/4' 



dm 



Am 



1 

if' 



d<tai/dm 
da™/dm 



(8) 



Detector acceptances and radiative corrections to the ini- 
tial particles are the same for non-resonant and J /ip con- 
tributions to ISR production of and cancel in the 
ratio. The total non-resonant cross section includes FSR 
contributions, which we parameterize in terms of K, the 
ratio of the visible non-resonant total and ISR-only (FSR 
switched ofi^) cross sections. Using simulated events, we 
determine K = 1.11 ± 0.01 (statistical error only) for our 
selection criteria. 

The result of the fit is shown in Fig. We find 
R = 18.94 ± 0.44 with per degree of freedom x^/'^ = 
122/144. This fitted value of R must be multiplied by 
1.002 to correct for non-resonant and resonant contribu- 
tions from e'^e~ — > e+e~7, 7r"'"7r~7, K^K~"f. The non- 
resonant cross section extracted from this measurement 



statistical error of K factor 0.9% 
systematic error of K factor 1.3% 
background uncertainty 0.5% 
simulation of J/ip line shape 1.4% 
interference effect 0.3% 



total 



2.2% 



TABLE I: The sources of systematic errors in R- K. 



is close to the value expected from simulation. Their ra- 
tio is 0.968 ± 0.016. The quoted uncertainty includes a 
0.4% statistical error, a 1% statistical error from simu- 
lation, and a 1.2% uncertainty in luminosity. We have 
not studied the systematic uncertainties on the efficiency 
for the non-resonant process in detail, as most of these 
cancel in R and, hence, do not affect the measurement of 
the J/tp parameters. The fitted value of Mj/.^ is shifted 
from that in the simulation by -(1.6±0.3) MeV/c^. The 
fitted value of is 3.4 ± 1.4 MeV/c^, corresponding to 
an overall mass resolution (sa 14.5 MeV/c^) 3% larger 
than that of the simulation. The background slope a 
corresponds to a 10% change of the non-resonant cross 
section in the mass range from 2.8 to 3.4 GeV/c^. The 
value of b is consistent with zero, in agreement with the 
Monte Carlo calculation. 

As seen in Eq. the ISR J/^j production cross sec- 
tion is proportional to the product of i?, determined from 
fitting the data, and K , determined from Monte Carlo 
simulations. The primary sources of systematic uncer- 
tainties for the product R - K are summarized in TableQ] 
Uncertainty in K is caused by different detection efficien- 
cies for the pure ISR process of J /ip production and the 
non-resonant e"*"e~ —t /i+/i^7 process to which both ISR 
and FSR amplitudes contribute. We estimate the uncer- 
tainty due to data-Monte Carlo differences by studying 
the stability oi R ■ K for different selection criteria. We 
vary the photon and muon angular selection criteria and 
the muon momentum requirement over a wide range of 
values. While the value of K varies from 1.08 to 1.19, 
the maximum deviation from our reference mean value 
R ■ K = 21.03 is only 1.3%. Although this variation 
might be a statistical fluctuation (at least in part), we 
treat it as a systematic uncertainty associated with the 
value of K . 

As described earlier, we use Monte Carlo simulations 
of specific ISR and other processes to estimate the level of 
non-resonant background to be less than 0.4%. We also 
use the data themselves to estimate this quantity. We 
compare the fit results for data selected with the standard 
selection criteria and for data selected with additional 
muon identification for one of the charged particles. This 
reduces pion (kaon) contamination by a factor of 9 (3). 
From the difference in i? • if, we estimate that the level 
of non-resonant background does not exceed 0.5%. 

The fit results do depend significantly on the model 
assumed for the J/ip line shape. The shape of the sig- 
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nal distribution varies with the selection of the maxi- 
mum allowed value of the ^ from the one-constraint fit. 
Requiring lower values tends to reject events with extra 
photons, thus reducing the fraction of events in the low 
mass tail of the J/V' peak. The fraction of J/'0 events 
with mass less than {Mjj^ — 0.1) GcV/c^ changes from 
2.4% for no cut, to 0.4% for < 5. Re-fitting data 
with different requirements on the value of does not 
change the result for R ■ K significantly. The maximum 
deviation oi R ■ K from our reference mean value, 1.4%, 
is taken as a systematic uncertainty. 

We also consider an additional contribution to the line- 
shape uncertainty by rc-fitting the data with a model that 
docs not include interference between the non-resonant 
and J/ip production amplitudes. The quality of this fit 
is good: x^ jv — 138/144. As the data do not distin- 
guish between the two models statistically, we take the 
difference in R, 0.3%, as the corresponding systematic 
uncertainty. The total systematic error for K -Rxs 2.2%, 
compared to the statistical error of 2.3%. 

The Born cross section for the process e+e~ 
J/ipl ~^ M^A*~7 can be evaluated from Eq. (jSJ. The 
non-resonant Born cross section in this formula is cal- 



culated to be 



d<T = 



4 MeV/c^ = 101.0 fb. Following 



the generally accepted practice [l^ of including the vac- 
uum polarization correction in the value of the electron 
width Fee, we multiply the pure Born cross section by 
1.042 ± 0.002. From R - K ^ 21.03 ± 0.49 ± 0.47 we cal- 



culate the cross section a 



j/i, 



2124±49±47fband the 



product of the J ftp parameters 



Fop ■ Bi, 



0.3301 ± 0.0077 ± 0.0073 keV. 



From the PDG values for B^e and B^^, which are dom- 
inated by those measured in 'tp{2S) — > J/i/'tt+tt^ decays 
by the BES Collaboration [^, we derive the electronic 
and total widths of the J/tjj meson, 

Fee = 5.61 ± 0.20 keV, F = 94.7 ± 4.4 kcV , 

using the correlated errors reported by BES. The statisti- 
cal and systematic uncertainties are combined in quadra- 
ture. Our results agree with the previous world aver- 
ages 0, Fee = 5.26 ± 0.37 keV and F = 87 ± 5 keV, but 
arc more precise. 
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